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Measurement of van-der-Waals interaction by atom trajectory imaging 
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We study the repulsive van der Waals interaction of cold rubidium 705'i/2 Rydberg atoms by 
analysis of time-delayed pair correlation functions. After excitation, Rydberg atoms are allowed to 
accelerate under the influence of the van der Waals force. Their positions are then measured using 
a single-atom imaging technique. From the average pair correlation function of the atom positions 
we obtain the initial atom-pair separation and the terminal velocity, which yield the van der Waals 
interaction coefficient Cq. The measured Ce value agrees well with calculations. The experimental 
method has been validated by simulations. The data hint at anisotropy in the overall expansion, 
caused by the shape of the excitation volume. Our measurement implies that the interacting entities 
are individual Rydberg atoms, not groups of atoms that coherently share a Rydberg excitation. 

PACS numbers: 32.80.Ee, 34.20.Cf 


The van der Waals interaction is important in the de¬ 
scription and control of interactions in few- and many- 
body dynamics studies. This interaction has been crit¬ 
ical in the observation of Rydberg excitation blockades 
and collective excitations HHa, Rydberg crystals EE], 
and Rydberg aggregates m m. Rydberg interactions 
have been used in quantum information processsing ler 
[12]. The van der Waals interaction between two Ryd¬ 
berg atoms has been measured using spectroscopic meth¬ 
ods m HU- Several measurements have been per¬ 
formed near surfaces to observe radiative Rydberg-level 
shifts caused by image charge interaction near metal sur¬ 
faces [m ng. The van der Waals interaction between 
excited cesium atoms and a dielectric surface has been 
measured using selective reflection spectroscopy m- 

Here, we develop a method to study the van der Waals 
interaction between Rydberg atoms using direct spatial 
imaging of their trajectories [THHUj. Pairs of 705'i/2 ru¬ 
bidium Rydberg atoms are prepared with a well-defined 
initial separation by detuning an excitation laser and uti¬ 
lizing the r~^ dependence of the van der Waals interac¬ 
tion After preparation, the atoms are subject to 

van der Waals forces (which are repulsive in this case). 
The effect of the forces is observed by tracking the inter¬ 
atomic distance between the Rydberg atoms, after they 
have been allowed to move for selected wait times (see 
Fig. [^. The atom trajectories and thereby the van der 
Waals interaction coefficient Cq are extracted from the 
pair correlation functions of the Rydberg atom positions. 

The experimental setup is shown in Fig. [^a). ®®Rb 
ground-state atoms are prepared in a magneto-optical 
trap (MOT) at a density of > 10^^ cm“^. The two- 
photon Rydberg excitation to 705'i/2 is driven by simul¬ 
taneous 780 nm and 480 nm laser pulses with a 5 /is 
duration and GHz red-detuning from the 5 P 3/2 inter¬ 
mediate state. Both beams propagate in the xy plane and 
are linearly polarized along z. The 780 nm beam has a 
Gaussian beam parameter wq of 0.75 mm and the 480 nm 
beam is focused to rco = 8 ym. The Rydberg atoms are 
ionized by applying a high voltage to a tip imaging probe 
(TIP) with a rounded tip of diameter 125 ym. Ions are 


accelerated by the TIP electric field towards a microchan- 
nel plate (MCP). Ion impacts result in blips produced 
by the MGP-phosphor detector assembly, revealing the 
Rydberg atom positions. In each experimental cycle we 
record one image, taken by a GGD camera; typically it 
contains several blips. For more detail see reference [21]. 

The excitation volume is ?^470 ym above the tip, 
which, in combination with the radial divergence of the 
TIP electric field and the MGP front voltage, results in a 
measured magnification of 155 with an uncertainty of 2%. 
To calibrate the magnification, we generate a Rydberg- 
atom distribution with a known spatial period by imaging 
the 480 nm beam through a double slit into the excita¬ 
tion region. The spatial period of the signal on the MGP 
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FIG. 1. (Color online) Geometry (a) and timing sequence 
(b) of the experiment. The wait time is the time difference 
between the end of the excitation and the onset of the field 
ionization pulse, (c) Van der Waals potential versus separa¬ 
tion between two Rydberg atoms. The separation increases 
due to the repulsive van der Waals force. We image the re¬ 
sulting trajectories by varying the wait time. In the insets 
in (c), circles represent two Rydberg atoms evolving under 
the influence of the van der Waals force. The (final) kinetic 
energy release, indicated by the blue velocity arrows, equals 
the van der Waals potential Wo at the initial separation ro. 
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FIG. 2. Pair correlation images from experiment (top row) 
and simulation (bottom row) for the selected wait times. The 
linear grayscale ranges from 0 (white) to 2 (black). Values 
of 1, <1, and >1 indicate no correlation, anticorrelation, and 
positive correlation, respectively. 


then yields the magnification. The resolution in the ob¬ 
ject plane is ^ 1 /im; it follows from the magnification 
and the point spread function of the detector. We can 
easily resolve the correlation length between the Rydberg 
atoms, which is on the order of 10 jam. The field of view 
typically contains 5 to 10 detected atoms. In the image 
processing sequence, we first extract ion impact positions 
on the MCP by using a peak detection algorithm [21]. 
Out of 10000 images taken in each dataset, we select the 
5000 images with the highest numbers of detected ions. 
From this subset, we calculate the average pair correla¬ 
tion image and normalize it such that at large distances 
it approaches the value of one. 

The excitation of the 705'i/2 Rydberg level is detuned 
by = 4 ± 2 MHz with respect to the two-photon reso¬ 
nance. The positive excitation detuning is compensated 
by the (repulsive) van der Waals interaction. We prefer¬ 
entially excite Rydberg atoms in pairs at a separation at 
which the van der Waals interaction is 26ui, = 8±4 MHz. 
We determine ro from pair-correlation images. After ex¬ 
citation, the Rydberg atoms are allowed to move within 
selected times before their positions are measured by ap¬ 
plying a field ionization pulse [Fig. [^b)]. The atom- 
position data are processed as explained above, yielding 
average pair correlation images. In Fig.|^ we present the 
pair correlation images for wait times 2.5 /is, 10 /is, 20 /is, 
30 /is, and 40 /is. The average interaction time of the 
Rydberg atoms is the wait time plus half the excitation 
pulse length (wait time + 2.5 /is). The pair correlation 
image at 2.5 /is exhibits strong correlation enhancement 
at a fairly well-defined radius. The initial correlation 
is critical for our trajectory experiment; it is sufficient 
to track Rydberg-pair trajectories out to approximately 
40 /is. With increasing wait time, the radius of enhanced 
pair correlation increases, reflecting an increase of inter¬ 
atomic separation due to the repulsive van der Waals 
interaction. At long wait times, the enhancement ring is 
blurred out due to initial thermal atom velocities in the 
MOT (temperature ^ 100 jiK). 

To determine the most probable separation rp(t) be¬ 
tween Rydberg-atom pairs at wait time t, we first com¬ 
pute the angular integrals I{r) of the experimental pair 


FIG. 3. (Golor online) (a) Angular integrals I(r) of the pair 
correlation images in Fig. at wait times: 2.5 /xs (red solid), 
10 /xs (orange dashed), 20 /xs (green dotted), 30 /xs (blue dot- 
dashed), and 40 /xs (pink short-dashed). The y axis is for the 
red-solid curve; for clarity, the other curves are shifted down 
in equidistant intervals of 0.3. (b) Interatomic separations 

between Rydberg-atom pairs as a function of interaction time 
(left axis, blue circles), obtained from the peak positions in 
the I(r) curves. The blue dashed line represents a linear fit 
at long wait times (20-40 /xs). The red solid line shows sim¬ 
ulation results obtained for = 3 MHz. The pink hollow 
squares show the visibility of the experimental pair correlation 
enhancement (right axis) as defined in text. 

correlation images in Fig. as shown in Fig. [^a). The 
values of r^{t) are obtained from local parabolic fits cen¬ 
tered approximately at the peak positions of the /(r) 
curves. We include 4 to 7 data points in the fit (depend¬ 
ing on the shapes of the curves). The resulting separa¬ 
tions rp(t) are shown in Fig. [^b). The visibility of the 
pair correlation enhancement, also shown in Fig. ib), 

IS (/max - /min)/(/max + /min), whcrC /^ax IS thc peak 

value of /(r) in the range > 10 /xm and /min is the min¬ 
imal value near 7 /xm. A cursory inspection of Fig. |^b) 
already shows that the trajectory of the Rydberg-atom 
pairs is characterized by an initial acceleration phase, 
during which the initial van der Waals potential energy, 
Wo, is converted into kinetic energy, and a later phase 
during which the atoms keep separating at a fixed veloc¬ 
ity. The drop in visibility is due to the thermally-induced 
blurring of the correlation ring at late times. 

To extract the van der Waals Cq coefficient, one may 
consider an isolated atom pair excited at an initial sepa¬ 
ration ro- The initial van der Waals energy is 

Wo = ^ = 2 ( 5 i / l . ( 1 ) 

ro 

From Eq. Ce = 2(5 z^l^o can, in principle, be obtained 
from spectroscopic measurement [14]. This method re¬ 
quires a well-defined ro, a narrow laser linewidth, and 
accurate knowledge of 6ui^. In our work, the relative un¬ 
certainty in Vq is about 12% (because the relative mag¬ 
nification uncertainty is 2%). In comparison, the uncer¬ 
tainty arising from Sui, =4 + 2 MHz is much larger. Tak¬ 
ing all uncertainties into account, Eq. [pleads to values of 
Cq ranging from 7x10“^^ Jm^ to 40x10“^^ Jm^. There¬ 
fore, Eq. only allows us to perform a crude order-of- 
magnitude estimate for Cq . The main weakness of Eq. 






3 


is that the Ce-values derived from it directly reflect our 
large relative uncertainty in 6 ui,. 

A better way for us to obtain Cq is to track the evolu¬ 
tion of the Rydberg-atom trajectories. Over the exper¬ 
imentally investigated wait times, the entire initial van 
der Waals energy Wq becomes converted into kinetic en¬ 
ergy. With the reduced mass of the atom pair, and 
the terminal relative velocity, Ut, it is 

= ^ = \livl ( 2 ) 

To obtain ro, we take the weighted average of the fit re¬ 
sults rp(t) at the earliest times used in the experiment 
[the points in Fig. [^b) at 250 ns and 2.5 /is]. The aver¬ 
aging is valid because during the first few microseconds 
after excitation the Rydberg atoms are frozen in place 
due to their inertia. The statistical weights are given by 
the inverse squares of the fitting uncertainties of at 
250 ns and at 2.5 /is. We obtain ro = 8.78 /im with a net 
fitting uncertainty of 0.02 /im. To determine Ut, we per¬ 
form a linear fit at long wait times (20-40 /is) and obtain 
Ut = 0.182 m/s with a fitting uncertainty of 0.008 m/s. It 
follows Cq = = (5.4T0.5) x 10“^^ Jm^. Including 

the 2% magnification uncertainty, the total uncertainty 
of Cq becomes 1.0 xl0“^^ Jm^. 

The final relative uncertainty of Cq follows from three 
statistically independent contributions: the magnifica¬ 
tion uncertainty, the fit uncertainty for Ut, and the fit 
uncertainty for vq. The respective powers at which these 
quantities enter into Cq are 8, 2, and 6. Factoring in 
these powers, the three quantities contribute respective 
independent uncertainties of 16%, 8.9%, and 1.4% to the 
relative uncertainty of Ce, leading to the total relative 
uncertainty of 18%. In our method, achieving a small 
magnification uncertainty is particularly important. 

A semi-classical 3D simulation of the dynamics of Ryd¬ 
berg atoms interacting by an isotropic van der Waals 
force has been performed to interpret the experiment. 
The simulation volume of (140 /im)^ exceeds the exper¬ 
imentally relevant volume by about a factor of two in 
each dimension. The 480 nm excitation beam propagates 
along the y direction. The excitation volume in the trans¬ 
verse directions {x and z) is limited by the size of the ex¬ 
citation beam {wq = 8 /im). The number of simulations 
is 5000, in keeping with the number of images analyzed in 
the experiment. Rydberg-atom positions and velocities 
are initialized as explained below, and the positions and 
velocities are then propagated using a Runge-Kutta inte¬ 
grator that includes all pair-wise interatomic forces. We 
use Cq = 5.7 X 10“^^Jm^ m in the initialization and 
integration procedures. To avoid edge effects, we crop 
the xy processing area in the images to (70 /im)^ before 
calculating the average pair-correlation images and the 
radial functions I{r). 

In each simulation, we start by drawing the number 
of Rydberg atoms to be used, nRyd, from a Poissonian 
distribution. Initial trial positions of the Rydberg atoms 
are then drawn from a probability distribution that is 


Gaussian along x and z, with a rcg of 8 /im, and uniform 
along y (in close analogy with the experiment). For an 
atom i at a trial position r^, the effective detuning ^eff (i*i) 
due to all other atoms j at positions Yj which are already 
excited into the Rydberg state is 

Since the detuning is substantial, the first pair of 
Rydberg atoms is simultaneously excited via off-resonant 
excitation [21] (because the intermediate state in which 
there is only one Rydberg atom present is off-resonant). 
Therefore the value of 8 in Eq. for the first Rydberg- 
atom pair is set to be twice the laser detuning dvi,. For 
the excitation of additional atoms we use 8 = 81 ^ 1 , to 
simulate the stepwise addition of those atoms, which 
can be a near-resonant process with other atoms already 
present |8] . Detunings due to the Doppler effect are about 
300 kHz and are neglected. The excitation probabil¬ 
ity Pex(^eff) is a Gaussian centered at (5eff = 0 with a 
FWHM of 4 MHz (given by the excitation bandwidth) 
and Pex(^^eff = 0) = 1. A new Rydberg atom at position 
Yi is created if the excitation probability is larger than a 
number randomly drawn between 0 and 1. This proce¬ 
dure is repeated with new trial positions until the desired 
number of riRyd Rydberg atoms has been reached. The 
initial center-of-mass velocities of the atoms are assigned 
using a Maxwell distribution at temperature 100 /iK. 

In the simulation, we record the atom positions at the 
same wait times as used in the experiment. The resulting 
pair correlation images are shown for the case 8 vi, = 
3 MHz in the bottom row in Fig.[^ Black pixels along the 
left and right edges of each pair correlation are an artifact 
due to the normalization used in the image processing. 

From the simulated pair correlation images we cal¬ 
culate the I{r) curves, as in the experiment, in order 
to obtain the most probable separations r^(t) between 
Rydberg-atom pairs at each interaction time [see solid 
curve in Fig. ib), which is for = 3 MHz]. The 
simulated and experimental results for r^{t) are in good 
agreement. In order to test how well our experimental 
procedure reproduces the Cq coefficient that underlies 
the atomic interactions, we evaluate the simulated results 
for rp(t) using the same method that we also employ for 
the experimental data. The results for rg, Vt and Cq 
extracted from the simulated r^{t) curves are shown in 
Table [T| for four choices of 8 vi,. These Cq values agree, 
within the uncertainties, with the value that has been 
entered as a fixed input into the simulation. This finding 
validates the experimentally used procedure; in particu¬ 
lar it is seen that the method is not very sensitive on 8 vi,. 
The experimental and calculated m Cq values are also 
included in Table [T] for reference. 

The experimental, simulated, and expected values for 
Cq in Table [l] are in reasonable agreement. Moreover, we 
observe that the entire simulated curve rp(t) for 8 ui^ = 
3 MHz matches the experimental result very well [see 
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TABLE L Comparison of calculated, experimental, and sim¬ 
ulated results for Ce, ro, and vt. 



(/xm) 

Vt 

(m/s) 

Cs 

(xlO"®®Jm®) 

Calculation 

Experiment 

8.8±0.2 

0.182±0.009 

5.7 

5.4±1.0 

Simulations 

(^z/L = 2 MHz 

8.85±0.06 

0.17±0.02 

4.9 ±1.0 

(5z/l = 3 MHz 

8.14±0.03 

0.233±0.005 

5.5T0.3 

(5z/l = 4 MHz 

7.76±0.05 

0.278±0.003 

6.0T0.3 

(5z/l = 5 MHz 

7.51±0.04 

0.28±0.01 

4.9T0.5 




FIG. 4. Pair correlation function at wait time 30 fis from 
(a) experiment and (b) simulation at = 3 MHz. The 
anisotropic expansion of the atom ensemble causes the radius 
of enhanced pair correlation along the transverse direction (x) 
to be larger than along the excitation beam direction (y). 


Figs. I andj^b)]. Overall, our findings lend credibility 
to our method of measuring Cq. 

We note that higher-order quadrupole-dipole and 
quadrupole-quadrupole interaction potentials, which 
scale as r~^ and are not important at the distances 
relevant in our work, at the current level of precision. We 
have verified this in a calculation of interaction potentials 
in which all terms up to the quadrupole-quadrupole in¬ 
teraction have been included [24] . 

A question of interest is whether the interatomic force 
is acting on “superatoms”, i.e. Rydberg excitations 
shared among a number of ground-state atoms, or on 
individual atoms. In our case, there are several tens of 
atoms within one blockade sphere. In our analysis, we 
have assumed that the effective mass of the interacting 
entities is half the rubidium atom mass, i.e. we have im¬ 
plied that the interacting entities are individual atoms 
and not super atoms. For our density and blockade ra¬ 
dius, the total mass of a super atom equals that of several 
tens of atoms. The agreement between the Cq coefficients 
in Table |l] demonstrates that the interacting entities are 
indeed individual atoms. The finding implies that, dur¬ 
ing the course of the van der Waals interaction, excita¬ 
tions within superatoms become projected onto individ¬ 
ual atoms, which are then ejected from the initial super¬ 
atom volumes. The phenomenon has been predicted in 
[25] for superatom clouds interacting via a dipole-dipole 
interaction. In our work, we arrive at a similar conclu¬ 
sion for van-der-Waals-interacting Rydberg atoms in a 
3D system. 

Close inspection of the pair correlation functions from 
the simulation presented in Fig. [^reveals anisotropic ex¬ 
pansion behavior: the radius of enhanced pair correlation 
is larger in x than in y direction. In several experimental 
data sets we see some evidence of anisotropic expansion, 
such as in the experimental result shown in Fig. Ba). 
The anisotropic expansion is due to the cylindrical shape 
of the excitation volume (not the interatomic interaction, 
which is isotropic). Since the blockade radius is close to 
the transverse size of the excitation region, at most two 
atoms can be excited side-by-side in x direction, lead¬ 
ing to free, unimpeded expansion along that direction. 
In contrast, more than two Rydberg atoms can be cre¬ 


ated along the y direction. Therefore, along y the ex¬ 
pansion is slowed down due to multiple-atom repulsion. 
In our experiment, the overall Rydberg-atom density is 
high enough to sometimes observe this effect. 

In summary, we have studied the trajectory of 
Rydberg-atom pairs interacting by repulsive, isotropic 
van der Waals interaction. We have extracted the Cq 
coefficient from the experiment and compared it with 
simulations and calculations, and have observed good 
agreement. The result implies that the interaction oc¬ 
curs between individual atoms, not superatoms. We have 
observed indications of an anisotropic effect in the ex¬ 
pansion, caused by the excitation geometry. Future work 
may involve atom-atom interactions that are anisotropic, 
such as the dipole-dipole interaction. 
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SUPPLEMENTAL MATERIAL 


We have used the model presented in detail in [i] 
to calculate the interaction potential between a pair of 
705 'i/ 2 rubidium Rydberg atoms with dipole-dipole inter¬ 
action terms only, as well as with dipole-dipole, dipole- 
quadrupole and quadrupole-quadrupole terms. The pro¬ 
jection of the electron angular momenta onto the inter- 
nuclear axis, M, can take the values 0 or ±1. For M=0 
there are two Rydberg-pair 705'i/2 — 705'i/2 levels, la¬ 
beled A and B, while for M=±l there is only one (see 
table). QQ means all interactions up to and including 
quadrupole-quadrupole terms are included. DD means 
only dipole-dipole terms are included. 

It is seen that the atom pairs have very similar shifts 
on all potentials, and that interaction terms beyond the 
dipole-dipole term are not important in the atom dis¬ 


tance range of interest (R>7 /im). The results also agree 
well with a perturbative level-shift calculation in [ii ]. 



M=0 Level A 

M=0 Level B 

M: 

=1 

R 

QQ 

DD 

QQ 

DD 

QQ 

DD 

/xm 

MHz 

MHz 

MHz 

MHz 

MHz 

MHz 

6 

18.215 

18.223 

19.046 

19.052 

18.424 

18.431 

7 

7.263 

7.268 

7.604 

7.609 

7.349 

7.353 

8 

3.267 

3.269 

3.421 

3.424 

3.305 

3.307 

9 

1.613 

1.614 

1.69 

1.691 

1.632 

1.633 

10 

0.858 

0.858 

0.898 

0.899 

0.868 

0.868 
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